Observation of a 2D Bose-gas: from thermal to quasi-condensate to superfluid 
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We present experimental results on a Bose gas in a quasi-2D geometry near the Berezinskii, 
Kosterlitz and Thouless (BKT) transition temperature. By measuring the density profile, in situ 
and after time of flight, and the coherence length, we identify different states of the gas. In particular, 
we observe that the gas develops a bimodal distribution without long range order. In this state, 
the gas presents a longer coherence length than the thermal cloud; it is quasi-condensed but is 
not superfluid. Experimental evidence indicates that we observe the superfluid transition (BKT 
transition) . 

PACS numbers: 03.75.Lm, 67.25.-d, 64.70.Tg 



One of the most fascinating aspects of a Bose gas 
in the degenerate regime is the role of dimensionality. 
A 2D interacting Bose gas is superfluid at low enough 
temperature [l|, 0]. However, by contrast to the 3D case, 
there is no long range coherence and the coherence de- 
cays as a power lawjl], At temperatures above the 
Berezinskii-Kostcrlitz-Thoulcss (BKT) transition tem- 
perature, the gas is not superfluid. Due to proliferation 
of free vortices, the quasi-condensate (QC) is fractured 
into small regions of nearly uniform phase, whose size, 
which corresponds to the typical length of the exponen- 
tial decay of the coherence, is larger than the thermal 
de Broglie wavelength (A = y / 2irh 2 /A/fcgT, where T is 
the temperature and M the atomic mass). For higher T, 
this size becomes smaller and approaches A, as the gas 
crosses over to the thermal phase. 

Experiments on 2D bosonic systems, such as two- 
dimensional 4 He films Q and trapped Bose gases 0, [B| , 
are able to show signatures of the BKT transition. Other 
systems, such as the superconducting transition in arrays 
of Josephson junctions Q and a two dimensional lattice 
of (3D) Bose- Einstein condensates 0], also exhibit a sim- 
ilar transition. Another interesting observation was in 
two dimensional spin polarized atomic hydrogen on liquid 
4 He Q where a reduction in three-body dipolar recom- 
bination (which is usually associated with condensation) 
was observed well above the BKT transition tempera- 
ture. This observation results from a reduction of density 
fluctuations, which corresponds to quasi-condensation [9j 

In this letter, we present evidence of transitions in a 
quasi-2D Bose gas from thermal (normal gas), to quasi- 
condensate without superfluidity, to superfluid quasi- 
condensate (BKT transition). We explicitly identify the 
theoretically expected non-superfluid quasi-condensate, 
a feature not clearly seen in other experiments on a 2D 
trapped Bose gases [3, We use an interferometric 
method to study the coherence of the gas down to dis- 
tances smaller than the thermal de Broglie wavelength. 
Our results can be understood using the local density 



approximation (LDA) on a model [111 ] of a homogeneous 
system. More recently, calculations for a trapped sys- 
tem have been carried out using classical-quantum field 
methods [l2[ and quantum Monte Carlo methods (l3j . 

The BKT transition occurs at a universal value of the 
superfluid density n s = 4/A 2 . However, the total density 
at the transition, ?ibkt, is not universal and depends 
on the strength of the interactions [l[. Interactions in a 
weakly interacting Bose gas trapped in a quasi-2D geom- 
etry (trap with tight confinement, at large frequency u>o, 
along one axis) are described by a dimensionless coupling 
constant g — \/8ira/lo [l4j|, where a is the 3D scattering 
length and l = yJh/MvQ the extent of the ground state 
along the tight direction. In the limit of weak interactions 
(ln(l/<?) < 1), tibktA 2 ps ln(C*/<?) [HGH. Monte-Carlo 
simulations that calculate the density as a function of the 
chemical potential give C ps 380 UJJ. We can use those 
simulations for our trapped system by applying a LDA 
(see Ref. [ll|). This is valid even near the BKT transition 
if the healing length defined at Tbkt, r c — 2X/y/2ng, is 
smaller than the typical size of the cloud (from the chem- 
ical potential given in [ll|, we find that the Thomas- 
Fermi radius of the cloud scales as R « ^/gRx, where 
R T = y / 2k B T/muj'j_ is the size of the thermal cloud, io± 
being the frequency of the trap in the 2D plane). This 
implies that to satisfy the LDA and thus have a BKT 
transition at a higher temperature than the true Bosc- 
Einstein condensation (BEC) transition (l7j . one must 
satisfy Hui± < gk B T. For our typical experimental pa- 
rameters, g ~ 0.02, T ~ 100 nK and uj ± /2tt ~ 20 Hz, so 
the LDA is reasonable. 



Applying the LDA to the calculation of Refs. [ll|, we 
get (Fig. [T]) the total density at the center of the trap as a 
function of the total number of atoms TV (normalized to 
the BEC critical number for 2D trapped, non-interacting 
atoms, iVcrit = tt 2 /6 (k B T/hLu ± ) 2 When TV ap- 

proaches iVcritj the central density increases rapidly and 
the calculation predicts the appearance of a bimodal dis- 
tribution (Fig. [l] inset). This increase coincides with 
a reduction in density fluctuations, which Refs. llj de- 
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FIG. 1: Calculated total density (Tot), quasi-condensate 
(QC) density and superfluid (SF) density at the center of the 
cloud in a trap as a function of the total number of atoms. 
Inset: density profile for different numbers of atoms, from 
AT/iVcrit ~ 0.9 (solid line) to N/Ncrit ~ 1.1 (dotted line). 

fine as the appearance of a quasi-condensate. The BKT 
phase transition occurs when the density at the center 
reaches ?xbkt- For T > Tbkt (n < iibkt), the 2D 
cloud will be broken up by the free vortices into multi- 
ple, phase independent, local condensates (non supcrfuid 
quasi-condensate). For T < Tbkt (n > %kt), there 
will be a superfluid (quasi-)condensate with phase fluc- 
tuations (due to phonons and bound vortex pairs). We 
do not expect to see a dramatic change of the in situ 
image at the phase transition, because it affects mainly 
the phase coherence and not the spatial density profile. 
However, in our experiment, phase coherence is revealed 
in time of flight (TOF) or by interferometry. By all these 
methods, we see the progression from thermal to quasi- 
condensate to superfluid. 

Our experiment uses sodium atoms confined in a sin- 
gle, horizontal, quasi-2D, optical dipole trap, similar to 
that described in [ljjj]. The infrared (1030 nm) trap laser 
beam, is spatially filtered through an optical fiber and 
collimatcd. The beam is then focused into a sheet of light 
using a cylindrical lens and projected onto the atomic 
cloud using a telescope. The sheet of light has a waist 
(1/e 2 radius) in the vertical direction of » 9 zxm, and 
a Rayleigh range of 440 fjxci. The waist in the horizon- 
tal direction perpendicular to the propagation axis of the 
laser is 800 /zm. However, due to a small distortion of the 
laser intensity profile in this direction, atoms see a trap- 
ping frequency higher than it would be for a Gaussian 
profile with that waist. This results in a trapped cloud 
of atoms that is approximately a circular disk with an 
aspect ratio of 50. For a typical value of the laser power 
(500 mW), the trapping frequencies are 1 kHz and 18 Hz 
in the vertical and horizontal directions, respectively. 

We initially confine sodium atoms in a magnetic trap 
and cool them using rf- induced evaporation [201 ] - Before 
reaching quantum degeneracy, we adiabatically transfer 
the atoms to the optical dipole trap. After 5 s of evapo- 
ration in the dipole trap, we reduce the depth of the trap 
(initially around 5 fiK) during 1 s to a final value rang- 
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FIG. 2: Density profile of the cloud after 5 ms TOF. The two 
images (with the same axis scale, but different color scale) 
are taken for two different numbers of atoms. The graph 
represent the cross section of the cloud along a line passing 
through the center of the cloud along the direction of the trap 
laser (arrow) . The dot-dashed (red) line is a fit to the sum of 
two Gaussians, while the solid (blue) line is the wider of the 
Gaussians. 

ing from 5 fiK to 2 fiK and then hold for 5 s to complete 
the evaporative cooling. We set the final temperature by 
choosing the final trap depth. We vary the phase space 
density keeping temperature constant (same trap depth) 
and controling the initial number of trapped atoms. For 
a trap with a vertical confinement of ujq/2tt = 1 kHz the 
temperature is around T = 100 nK, which corresponds 
to ksT ~ 2fkuo, so the thermal motion is not completely 
frozen out along the tight direction. However, the quasi- 
condensed atoms are well confined in 2D since the mean 
field never exceeds 1 kHz, and therefore we expect to see 
2D physics. 

Figure [5] shows two absorption images of the cloud, af- 
ter 5 ms of TOF, for the same temperature but different 
number of atoms. Both pictures show bimodality and 
density fluctuations. These fluctuations, which appear 
only after TOF, are manifestations of the in situ phase 
fluctuations of the quasi-condensate. In spite of the sim- 
ilarities (bimodality and density fluctuations), there are 
significant differences between the two pictures. The one 
on the left with fewer atoms is clearly broken up, while 
the one on the right is more compact and the narrow part 
of the fit (dot-dashed/red curve) is not as wide. 

In figure [3] we plot the width of the narrow part as a 
function of the peak density for two different tempera- 
tures (trap depths). For both temperatures, the width 
initially decreases with increasing density (presumably 
due to a decrease of free vortices) and then, beyond a 
certain point, increases slowly (presumably due to repul- 
sive interactions between atoms). The central density at 
which the width is a minimum depends on the tempera- 
ture of the gas: we identify this as the point when atoms 
at the center of the trap cross the BKT transition. In 
the inset of Fig. [3J we plot this critical density (2l| as 
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FIG. 3: Width of the narrow part of the cloud, obtained by 
fitting 5 ms TOF images with two Gaussians, as a function 
of the fitted peak 2D density, for two different temperatures. 
The lines are linear fits to two portions of the data, and are 
used to determine a transition point. Inset: 2D "in situ" 
peak density [5lJ at the transition point as a function of the 
@hyphenationtem-pera-ture. The dashed line is the theoret- 
ical value tibktA 2 = ln(380/g) and the solid line is the same 
value corrected using a 3D model [2^]. We add a ±15% uncer- 
tainty [U primarily due to the disagreement of the different 
methods used for calibration (for the temperature: TOF and 
interferometric measurement of A; for the density: theoreti- 
cal absorbtion cross section and in situ measurement of the 
Thomas- Fermi radius). 
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FIG. 4: Cross section of the density profile of the cloud after 
10 ms of TOF, fitted with the sum of three Gaussian. Left: 
for n < riBKT (T > Tbkt); right: for n > jibkt (T < Tbkt). 
The data are taken with the same trap depth as Fig [2] 



a function of the tem per ature. The dashed line is the 
theoretical prediction ll| tibktA 2 = ln(380/<?) and the 
continious line is this theoretical prediction corrected by 
a factor taking into acount thermal excitation of modes 
in the tight direction, as done in Ref. (22|. The position 
of the point of minimum width agrees with the theoret- 
ical prediction, supporting its identification as the BKT 
transition in the central region of the trap. 

With this interpretation, the picture on the left of 
Fig. [2] is taken for T > Tbkt and contains a thermal 
component (solid/blue curve) and a non superfluid quasi- 
condensate (dot-dashed/red curves) whereas the picture 
on the right (T < Tbkt) contains the thermal part and a 
QC part with both a superfluid and non-supefluid com- 
ponent, which are not separated after 5 ms TOF. How- 
ever, after 10 ms of TOF, the three components can 
be separated. For T just above Tbkt, there is a bi- 
modal distribution (Fig. 0] left): the thermal component 
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FIG. 5: Left: Schematic of the atom interferometer used to 
measure the spatial coherence of the 2D atomic cloud. The 
us and ks denote the frequencies and corresponding wavevec- 
tors, respectively, of the Raman laser beams. Right: Im- 
ages of atoms in the F = 2 state after the interferometer 
sequence for two different delays between Raman pulses (up- 
per: r = 7.5 us, lower: r = 17.5 /is). For the longer delay, 
the interference of the thermal cloud is washed out and there 
are only fringes from the condensate. 



(solid/blue) and quasi-condensate (dot-dashed/rcd). For 
T just below Tbkt (Fig- HI right), we see a trimodal dis- 
tribution: the thermal component and the non-supcrfluid 
QC (which are very similar to the case where T is just 
above Tbkt) and the superfluid part (very narrow peak, 
dashed/orange). This transition occurs when the central 
density reaches «4.5 atoms//xm 2 , but with the approx- 
imate correction factor to the in situ value gives «9.4 
atoms//xm 2 , which is in reasonable agreement with the 
theoretical prediction and the 5 ms TOF analysis. 

To measure the coherence of our gas, we have devel- 
oped a Ramsey-like method. This allows us to measure 
coherence properties on a length scale below the ther- 
mal dc Broglic wavelength. We use two, 7r/2, two-photon 
Raman pulses, with counterpropagating and nearly coun- 
terpropagating laser beams, to transfer atoms from the 
F = 1 to the F — 2 internal state, also transferring two 
photon recoil momenta. We interfere, in situ, the two 
F = 2 copies of the cloud, which have a velocity differ- 
ent by KK /m and arc spatially shifted by R (see Fig. [5] 
and below). Specifically, if "0o(r) is the quantum field 
operator, with no(r) = ^g(r)^o( r )^ the initial atomic 

density (in F = 1), then the measured density of atoms 
in F = 2 after the two Raman pulses is 



n(r) = - (2no(r) + ($(r)^(r- 



R)e 



h.c. 



(i) 

where K = ki — k2 — k' x + k 2 (notation defined in Fig. [5]) , 
R = h(ki — k2)r/m is the displacement of the cloud 
during the time between pulses, r, and 4>o is the (un- 
controlled) phase between the two Raman pulses. By 
measuring the average spatial contrast of the interfer- 
ence fringes in an image, averaged over several images, 
we obtain the normalized correlation function g^{R) = 

$(r)$,(r-R))/no(r). 

Figure [5] shows the contrast of the fringes as a function 
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FIG. 6: Contrast of the interference fringes as a function of 
the separation of the two clouds in the interferometer. Each 
curve corresponds to an average over many data sets where 
the maximum density of atoms falls within a chosen range. In- 
set: Corresponding density profiles of the atomic cloud. Lines 
(solid, dotted, dashed, dot-dashed) correspond respectively to 
density ranges of 2-4, 5-8, 10-12, 13-19 atoms//im 2 . 



of R, for different densities, but the same temperature, 
T w 100 nK. This contrast is calculated by measuring 
the relative intensity of the peak of the Fourier compo- 
nent of the image at the spatial frequency of the fringes. 
Since we cannot precisely control the number of atoms, 
we average only the data that falls within a chosen range 
of peak densities. For each curve, we see an initial drop 
of the contrast on a length scale of the order of 1 /im, 
due to the relatively short coherence length of the ther- 
mal component (the thermal de Broglie wavelength is 
A = 1 /im for T = 100 nK). For atomic clouds with den- 
sities between 5 and 8 atoms/^m 2 (dotted line), where 
we clearly see a bimodal distribution, the coherence ex- 
tends well beyond A, but decreasing to zero by 10 /im, 
a distance much shorter than the spatial width of the 
narrow part of the bimodal distribution (see inset). This 
behavior is as expected for this quasi-condensate region 
for T just above Tbkt where, due to local condensation, 
there is coherence for distances greater than the thermal 
de Broglie wavelength, but because of the free vortices, 
the coherence decreases exponentially on a length scale 
2A/ A 



llj . For higher [1] 



of the order of r c = 2X/y/2^g ~ 6 /im 
densities (T < Tbkt), the coherence also decays slowly 
but has a non zero value even at 20 //m. For this phase, 
one expects a power law decay of the coherence mainly 
due to long wavelength phonons. The typical distance at 
which the coherence decreases by a factor of 2 is given 
by I = 2 1 / Q /v^ Q, where a is the coefficient of the 
power law decay whose value is 1/4 at the BKT transi- 
tion. We calculate I m 40 /im for our typical parameters 
at the BKT transition. 

With our interferometric method, an isolated vortex 
would be characterized by a double fork structure in 



a distance which is comparable to the healing length 
(r c ~ 6 /im). However, in the non-superfluid region 
(T > Tbkt) j the fringes disappear at this separation. 
Therefore, with our fringe spacing (constrained by our 
optical resolution), it is not possible to see such a direct 
signature of the free vortices that are supposed to prolif- 
erate in this regime. This should also holds true in 0], 
as suggested by the results of 

In this experiment, we have observed the predicted 
intermediate regime, between the thermal and super- 
fluid phase, characterized by a bimodal density distri- 
bution and a coherence length longer than the thermal 
de Broglie wavelength, but smaller than the predicted 
length scale in the superfluid region. At higher densi- 
ties, we see a phase with a much longer coherence length, 
which is consistent with the BKT superfluid phase. For 
sufficiently long TOF we clearly see a trimodal distribu- 
tion showing the presence of three components. A pos- 
sible explanation for previous experiments 0, Q having 
missed seeing a trimodal distribution (they associate the 
BKT transition with the appearance of bimodality) is 
that, they don't see this intermediate quasi-condensed 
regime. This may be because their imaging integrates 
over one dimension, reducing the signal of the narrow 
part of the distribution relative to that of the thermal 
cloud. Finally, our expermental setup, which allows us to 
sec BKT related physics in a single "pancake" of atoms, 
opens new oportunitics, such as studying density fluctu- 
ations at the transition. Future experiments could use a 
Fcshbach resonance to tune the coupling constant in our 
all optical trap, and explore different interaction regimes. 

This work was suported by the Office of Naval Re- 
search. We thank J. Dalibard, T. Simula, M. Davis and 
B. Blakic for helpful discussions. 



the interference pattern [24(. To see such a structure 



one should separate the two interfering copies of the gas 
by a distance larger than the fringe spacing (w 7 /im), 
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